Wringing out DNA 



O 
O 

tin 



Timothee Lionnet^Q Sylvain Joubaud^, Richard Lavery^, David Bensimon^, and Vincent Croquette^ 
^ Laboratoire de Physique Statistique, CNRS UMR 8550, 
Ecole Normale Superieure, 24 rue Lhomond 75005 Paris, France 
^ Laboratoire de Physique CNRS UMR 5672, Ecole Normale Superieure de Lyon, 46 allee d Italic 69364 Lyon Cedex 07, France 

^ Laboratoire de Biochimie Theorique, CNRS UPR 9080, 
Institut de Biologic Physico-Chimique, 13 rue Pierre et Marie Curie, Paris 75005, France 

(Dated: February 9, 2008) 

The chiral nature of DNA plays a crucial role in cellular processes. Here we use magnetic tweezers 
to explore one of the signatures of this chirality, the coupling between stretch and twist deformations. 
We show that the extention of a stretched DNA molecule increases linearly by 0.42 nm per excess 
turn applied to the double helix. This result contradicts the intuition that DNA should lengthen 
as it is unwound and get shorter with overwinding. We then present numerical results of energy 
minimizations of torsionally restrained DNA that display a behaviour similar to the experimental 
data and shed light on the molecular details of this surprising effect. 

PACS numbers: 87.14. Gg, 87.15.Aa, 87.15.La 
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The helical structure of double-stranded DNA (ds- 
DNA) results in very specific mechanical properties. 
These affect the function of nucleic acids, as they de- 
termine the accessibility of the genetic material to the 
proteins that process it. Micromanipulation techniques 
developed over the past decade have made it possible 
to manipulate a single DNA molecule, and thus to ex- 
tensively study its mechanical response Q, In the 
so-called entropic regime {F < 10 pN), stretched dsDNA 
behaves as a uniform, semi-fiexible rod whose extension 
results from a balance between the entropy of bending 
fiuctuations and the work performed by the stretching 
force 0, 0, 0, Q . In the very low force regime (typically 
F < 0.4 pN), twisted DNA behaves as an isotropic rod. 
Its torsional energy increases up to a threshold where the 
molecule buckles to form superhelical loops called plec- 
tonemes fsl, H llO, 111 IH 111. At higher forces (and 
consequently higher torques), DNA chirality affects its 
response to twist. The molecule undergoes structural 
transitions away from its native conformation in solu- 
tion, 5-DNA. When the torque is negative, i.e. opposite 
to the handedness of the double helix, DNA denatures: 
its unwound strands locally separate. When the torque is 
positive, the molecule adopts a highly overwound state, 

p-DNA jllillii. 

As a result of its helical structure, a DNA molecule is 
expected to display a coupling between its extension L 
and its degree of twist Tw. In other words, one expects 
the extension of DNA in its native structure to be an 
asymmetric function of twist. The response of DNA to 
both twist and stretch has been studied in the framework 
of continuum elasticity fTEO,E3- In that approach, the 
energy per unit length E/Lq is a quadratic function of 
its relative change in extension e = L/Lq — 1 and twist: 



a — Tw/Twq — 1 = n/Two, (where Lq and Twq are the 
contour length and natural twist of the double helix and 
n is the number of turns added to a twisted DNA): 



Be^ + 2De(j) - Fe 
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where p = Lq/Tw = 3.6 nm is the DNA pitch and C ^ 
100 nm, B ^ 78 nm and D are respectively the DNA 
torsional, stretch and twist/stretch moduli Q, E3- By 
minimizing E with respect to e at fixed <j, one finds that 
the extension should vary linearly with the twist: e — 
ecr=o = —Da/B, equivalent to AL = —pDn/B. Values of 
D in the 12-20 nm range have been previously extracted 
from single molecule results [16, 17J. This corresponds to 
a decrease in length of the DNA molecule of 0.6-0.9 nm 
per added turn. However, these estimates must be taken 
with care as they were extracted from data pertaining to 
DNA adopting different structures, not necessarily the 
native 5-DNA form. 

In the absence of a reliable twist-stretch estimate, elas- 
ticity theory needs to be backed by a complementary, 
atomic-scale approach, in order to predict the value of 
the phenomenological parameter D. Geometry alone pro- 
vides a first "naive" estimate, if we assume that the dou- 
ble helix radius R = 1 nm and the arc length of the 
backbones Ni^p^/R^O'^ + z'^ are fixed (A^^p being the total 
number of base pairs, z = L/N^p the rise per base pair 
-the difference between adjacent base pairs along the he- 
lical axis- and 6 = 27rTw/N})p the twist per base pair 
-the angle formed by adjacent base pairs around the he- 
lix axis- see FigGJA). The length of the molecule should 
then increase as it is unwound: dL/dn = —Att^R'^/p ^ 
- 11 nm/turn (or D = 240 nm). A finer approach based 
on Monte Carlo simulations suggests a similar estimate 
(2nl |. Elastic constants calculated from Molecular Dy- 
namics (MD) simulations yield D values ranging from 4 
to 18 mn for different base sequences, which implies an 
average value of dL/dn = -0.5 nm/turn [2l|. All these 
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FIG. 1: A: Two adjacent base pairs are schematically drawn 
as two pairs of plates, connected by the two backbones (thick 
black lines). Notice the negative value of the shift in B- DNA, 
i.e. the displacement of the base pair from the helical axis 
towards the minor groove (marked by 'm'). 'M' marks the 
major groove. B: Principle of the magnetic tweezers (not to 
scale). 



estimates, in spite of important differences in their mag- 
nitude, agree on the twist-stretch coupling sign: DNA 
should lengthen as it is unwound, as suggested by the 
simple picture of wringing out a floorcloth. 

We now describe magnetic tweezer experiments on sin- 
gle DNA molecules that allow measurement of the value 
of the twist-stretch coupling of 5-DNA. While its sign 
is opposite to the "naive" expectations previously men- 
tioned, it is however in good agreement with atomic-scale 
numerical modeling of the response of DNA to twist. The 
modeling results further suggest that response of 5-DNA 
to torsion illustrates the same mechanical coupling as 
seen in the transition between 5-DNA and the form of 
DNA observed under low hydration conditions, 74-DNA 

Magnetic tweezer experiments - Magnetic tweezers al- 
low to monitor the change in extension of a single dsDNA 
molecule as it is twisted ^23*1 . In a custom-built flow 
chamber, we tether a single DNA molecule between a 
glass surface and a superparamagnetic bead (1.4 /im ra- 
dius, Dynal). A pair of permanent magnets placed above 
the sample generate a constant, vertical stretching force 
on the bead and lock its rotational motion around the 
vertical axis. The exerted force F and the rotation n of 
the bead can be set by respectively translating and rotat- 
ing the magnets (FigQ^). The bead is imaged at video 
rate (60 Hz) through a lOOx oil immersion objective, with 
a typical resolution of 2 nm. The force F exerted on the 
bead is computed from measurements of {Sx'^) and L us- 
ing the equipartition theorem: F = ksTL/ {Sx'^) [23\. 

The curves L{n) are obtained by rotating the end of a 
DNA molecule under a constant force F and averaging 
the measured extension L for each value of n over typi- 
cally 128 points. At high enough F and for small n val- 
ues (typically -0.01 < a < 0.02), rotation translates into 
a change of the DNA twist. This approximation, which 
neglects relaxation of torsional constraint through bend- 
ing fluctuations ("writhe"), is correct (when F ^ 2 pN) 
within a 10 % margin [llj. However, for larger \a\ values, 
buckling and/or structural transitions invalidate this re- 
lationship. This is illustrated in Figl^lA, where dsDNA 
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FIG. 2: A: Extension L versus rotation n curve of DNA 
(pLacDl, 7.4 kbp; 20 mM Tris pH 7.9, 25 mM NaCl, 0.625 
mM MgCl2, 0.1 % w/v BSA; F = 7 pN). In the 5-DNA sta- 
bility region (black overlay), the curve displays a linear be- 
haviour, with a slope of 0.28 nm/turn. B: Effect of the buffer 
conditions, exerted force and DNA sequence content on the 
twist-stretch coupling. Since the slope value does not display 
any significant variation, we extract its average value (S.D.) 
in nm/turn from each set of data: solid squares: 0.52 (0.13); 
black line: 0.39 (0.05); open circles: 0.70 (0.24); solid gray 
circles: 0.43 (0.21). 



is stable in its B- form for small torsional strains (black 
overlay, -0.01 < a < 0.02). Storing twist energy through 
small deformations around the relaxed B- conformation 
becomes unfavorable under higher torsional strain, and 
dsDNA undergoes a transition to supercoiled P-DNA (sc- 
P) in the positive supercoiling region (a > 0.02), whereas 
strong negative supercoiling conditions induce local de- 
naturation of the molecule (a < -0.01) ^,i2^J. The twist- 
stretch coupling measures the change in extension for 
small variations in twist. It is thus deduced from the 
slope of L{n) at n -0: ^)n=o = 0.28 ± 0.04 nm/turn 
(Fig|2lA) equivalent to a twist-stretch modulus D = -6.1 
± 0.9 nm. Notice that as B-DNA is overwound its exten- 
sion increases, in contradiction with one's naive intuition. 

The twist-stretch coupling is neither force- nor [Mg^+]- 
dependent (within the 2-20 pN and 0-10 mM range ex- 
plored (Fig|2}. We did not observe any sequence depen- 
dence, as two unrelated molecules with different sequence 
contents (pLacDl: 7.4 kbp, 43 % GC; pSA850: 3.6 kbp, 
48 % GC) yield similar values. By averaging all the val- 
ues obtained in the different conditions, we obtain a slope 
of 0.42 ± 0.2 nm/turn, and thus D = -9.1 ± 4 nm. 
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Molecular modeling calculations - In the absence of a 
satisfactory theoretical description, we decided to inves- 
tigate the atomic details of dsDNA response to twist in 
order to understand our experimental results. This has 
been done by performing energy minimization of DNA 
under conditions mimicking our experiments. 

The modeling calculations were carried out with 
JUMNA [25] using the AMBER parm 98 force field 
DNA was modeled as a helically symmetric polymer with 
a mononucleotide or dinucleotide symmetry repeat ac- 
cording to the nature of the base sequence. Solvent and 
counterion effects were modeled using a sigmoi'dal dis- 
tance dependent dielectric function and reduced phos- 
phate charges (-0.5 e). Calculations made with a gener- 
alized Born solvent model [27] gave very similar results. 
Sequence effects were investigated by making calculations 
for all dinucleotide repeating sequences, (AA)^, (GG)n, 
(AT)n, (CG)n, (AC)n and (AG)n, and for an "average" 
sequence composed of equal contributions from each of 
the four standard base pairs, obtained using the multi- 
copy approach ADAPT [28|. For each sequence, the twist 
per base pair was fixed at values ranging from 32.5° to 
40.5° (at intervals of 0.25°). The DNA structure was 
then energy minimized and its helical conformation was 
analyzed. The effects of tension were studied by applying 
equal and opposite forces to both ends of each strand of 
the double helix [29]. 

In the absence of an exerted force, the base pair exten- 
sion, or rise, z is a linear function of the twist per base 
pair in the positive supercoiling region, i.e. when the 
twist is greater than its relaxed state value Oq ^ 34.5° 
(FigEK). All tested sequences have comparable slopes, 
typically 0.9 nm/turn (thus D = -20 nm). In the nega- 
tive supercoiling region, the curves display different be- 
haviours that cancel out to produce a fiat line in the 
case of the "average" sequence. Exerting a force on the 
molecule slightly decreases the slope of the rise versus 
twist curve. At 6 pN (18 pN), the rise increases by 0.68 
nm/turn (0.58 nm/turn), equivalent to = -15 nm (-13 
nm). 

Other basepair parameters, such as inclination, roll, 
etc I30] were monitored as the twist is varied. As shown 
in FiglS^, increasing rise is tightly coupled to more neg- 
ative inclination {i.e. a counterclockwise rotation of the 
base pairs viewed from the minor groove of the double 
helix). The absence of significant rise variation observed 
for energy minimizations of a DNA molecule with locked 
base pair inclination further illustrates this tight coupling 
(data not shown). Other structural parameters display 
a significant change in this range of twist values: the 
shift becomes less negative with increasing twist, with a 
slope of 0.38 A/degree (shift, or, more formally, the X 
displacement is the displacement of the base pairs from 
the helical axis perpendicular to their long axis, and is 
positive for displacements towards the major groove); the 
sugar phase angle [22]) increases by 2.5 degree per degree 
of twist; finally, the double helix diameter decreases by 
-0.44 A/de erree. 



FIG. 3: A: Twist-stretch coupling extracted from energy min- 
imizations (F = 0), for different dinucleotide sequences. All 
curves display a linear behaviour around the relaxed state. 
Fitting each curve in its linear region to a line yields the value 
of the slope in nm/turn: 0.9 (average), 1.15 (GC), 0.97 (AT), 
1.33 (AC) and 0.76 (GA). a = refers to the relaxed twist 
per base pair for the average sequence (34.5°). B: variation 
of the rise (dashed) and inclination (solid) as functions of the 
induced twist for the GA sequence, a = refers to the relaxed 
twist per base pair for the GA sequence (36.75°). 



Concluding remarks - The magnetic tweezer experi- 
ments presented in this work provide clear evidence of a 
linear twist-stretch coupling in 5-DNA, in the regime of 
small, physiological torsional strains {\a\ < 0.05). The 
sign and magnitude of this effect differ from naive theo- 
retical predictions. However, they are in good agreement 
with the present molecular modeling calculations of the 
response of 5-DNA to limited changes in the twist. 

A few details remain to be adressed. Calculations pre- 
dict that D should decrease with increasing F. Exper- 
imental data lack the necessary resolution at low force 
to observe such an effect. In addition, the range of su- 
percoiling values in which this effect is observed differs 
between experiments and calculations. Energy minimiza- 
tions include neither the possible structural transitions to 
different DNA states (denatured or P) nor sense of buck- 
ling transition. This explains why a larger a range can 
be studied than in the tweezer experiments. However, it 
is not clear why linear, positive twist-stretch coupling is 
only observed for positive supercoiling in the calculations, 
whereas experimentally, such a behaviour is observed for 
both positive and negative supercoiling. 
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FIG. 4: Axial and lateral views of the structures obtained 
by energy minimizations (14 base pairs; "average" sequence). 
Only the phosphodiester backbones and base pairs are repre- 
sented. Left: relaxed 5-DNA state, 6* = 34.5°. Right: 5-DNA 
under positive supercoiling (Q = 40.5°) displays a smaller di- 
ameter, an increased length, a negative base pair inclination, 
and a less negative shift value. 



present results to MD simulations (2lL l3ll|. 

The present numerical results yield evidence for the 
molecular mechanism of the observed effect: the decrease 
in rise as the molecule is unwound is coupled to an in- 
clination of the basepairs towards positive values, an in- 
crease in the diameter, more negative shift values, and a 
decrease of the sugar phase angle (FigE)). All these vari- 
ations, in spite of a much smaller amplitude, have the 
same sign as the variations observed in the B- to 74-DNA 
transition [22j. This suggests that the twist-stretch cou- 
pling involves the same helical deformation mechanism 
as the transition to 74-DNA. This behaviour might be 
important as part of the mechanism of DNA deforma- 
tion which is exploited by proteins in detecting specific 
sequences ^3^ . 



The difference between our results and prior D esti- 
mates based on single-molecule experiments is most likely 
due to the presence of extreme structural transitions in 
the data previously exploited [16, 17J. This suggests 
that structural transitions exploit very different atomic 
mechanisms than fluctuations within the 5-DNA stabil- 
ity region. Work remains to be done to clearly relate the 
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